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State-of-the-art ab initio techniques have been applied to compute the potential energy curves for the elec- 
tronic states in the A-'-E^, c 3 II u , and a 3 S+ manifold of the strontium dimer, the spin-orbit and nonadiabatic 
coupling matrix elements between the states in the manifold, and the electric transition dipole moment from 
the ground X*E+ to the nonrelativistic and relativistic states in the A+c+a manifold. The potential energy 
curves and transition moments were obtained with the linear response (equation of motion) coupled cluster 
method limited to single, double, and linear triple excitations for the potentials and limited to single and 
double excitations for the transition moments. The spin-orbit and nonadiabatic coupling matrix elements 
were computed with the multireference configuration interaction method limited to single and double excita- 
tions. Our results for the nonrelativistic and relativistic (spin-orbit coupled) potentials deviate substantially 
from recent ab initio calculations. The potential energy curve for the spectroscopically active (1)0+ state 
is in quantitative agreement with the empirical potential fitted to high-resolution Fourier transform spectra 
[A. Stein, H. Knockel, and E. Tiemann, Eur. Phys. J. D 64, 227 (2011)]. The computed ab initio points 
were fitted to physically sound analytical expressions, and used in converged coupled channel calculations 
of the rovibrational energy levels in the A+c+a manifold and line strengths for the A 1 !^ <(— X 1 E+ tran- 
sitions. Positions and lifetimes of quasi-bound Feshbach resonances lying above the x So + 3 Pi dissociation 
limit were also obtained. Our results reproduce (semi)quantitatively the experimental data observed thus far. 
Predictions for on-going and future experiments are also reported. 



I. INTRODUCTION 



In recent years the strontium diatomic molecule, Sr2, 
has attracted the interest of theoreticians and experi- 
mentalists. Similarly to the calcium dimer, Sr2 in its 
ground X X E+ state does not form a chemical bond. 
Indeed, the binding energy of 1081.6 cm" 1 - is typical 
for weakly bound complexes rather than for chemically 
bound molecules. However, excited states of Sr2 are 
strongly bound, and have been observed in many exper- 
iments. As a matter of fact, the strontium molecule was 
the subject of numerous high-resolution spectroscopic 
studies in the gas phased 3 - and in rare gas matrices^. 
The dissociation energy of the ground state was first esti- 
mated from the RKR inversion of the spectroscopic data 
for the I^Ej «— X X E+ transitions 3 . Recently, a more 
elaborate study of Sr2 was reported^, in which measure- 
ments of the l^Ej «— X^-El" transitions covering large 
parts of the ground state well were recorded. The spec- 
trum corresponding to the B 1 E+ <— X X E+ transitions 
could easily be assigned using standard spectroscopic 
techniques because the B state dissociating into X S + 1 P 



atoms is relatively well isolated, thus not perturbed by 
any other electronic state. While the same is true for the 
A' U u state reported in Refi£, this is not the case for the 
A-'-E* state dissociating into 1 S + *D atoms. Here the 
potential energy curve of the A state crosses the curve of 
the c 3 n M state dissociating into *S + 3 P states, and the 
corresponding spectrum cannot easily be assigned. Ex- 
perimental investigation of the A state is limited to Rcf. 6 . 
Tiemann and collaborators measured the spectrum cor- 
responding to the A^^E,^ «— X-'-E* transition by high- 
resolution Fourier transform spectroscopy, and showed 
that the rovibrational levels of the A state get strongly 
perturbed by the c 3 H u state 6 . Unfortunately, the amount 
of the measured data was not sufficient to apply a de- 
perturbation procedure that could be used with trust to 
determine the spectroscopic constants of the A*E+ and 
c 3 n M states. In particular, the key information on the 
spin-orbit coupling between the A and c states could not 
be determined from the analysis of the spectra, and only 
an effective potential was obtained. 

Most of the ab initio calculations on the Sr2 molecule 
reported in the literature thus far are concerned with 



the ground state potential energy curved and the van 
der Waals constants governing the long-range behavior 
of the ground state potential 9 - - — . To the best of our 
knowledge only three theoretical papers considered the 
excited states of the strontium dimer— ~—. However, in 
view of the recent findings of Tiemann and collaborators 6 
the quality of these data is questionable. 

It should be stressed that alkaline-earth atoms and 
molecules are not only interesting for conventional spec- 
troscopy, but are also intensely investigated in experi- 
ments at ultralow temperatures. Closed-shell atoms such 
as alkaline-earth metal atoms are much more challeng- 
ing to cool and trap than open-shell atoms like the al- 
kali atoms. They do not have magnetic moments in 
the ground state that would enable magnetic trapping. 
Moreover, the short lifetime of the first excited 1 Pi state 
implies rather high Doppler temperatures, requiring a 
dual-stage cooling with the second stage operating near 
the 3 Pi intercombination line. Despite these challenges, 
cooling of calcium, strontium, and ytterbium atoms to 
micro-Kelvin temperatures has been realized, and Bosc- 
Einstein condensates of 40 Caii, 84 Sr.±&±2, 86 Sr^, 88 Sj21, 
170 Yb 22 , and 174 Yb 2 ^ have been obtained. 

On the other hand, the closed-shell structure of the 
alkaline-earth metal atoms leads to very simple molecu- 
lar potentials with low radiative losses and weak coupling 
to the environment. This opens new areas of possible ap- 
plications, such as manipulation of the scattering prop- 
erties with low-loss optical Feshbach resonances 24 , high- 
resolution photoassociation spectroscopy at the inter- 
combination lin o 25 ' 26 , precision measurements to test for 
a time variation of the proton-to-electron mass ratio 27 ' 28 
and of the fine structure constant^ 9 -, quantum compu- 
tation with trapped polar molecules 3 ^, and ultracold 
chemistry 31 . 

Of particular interest for the present work are the ex- 
perimental investigations of the photoassociation spectra 
near the intercombination lin o 25 ' 26 ' 32 and proposed pre- 
cision measurements of fundamental constant s 27 ' 28 . In 
particular the latter require a precise route for the pro- 
duction of ultracold molecules in predefined rovibrational 
states. This, in turn, requires an accurate knowledge 
of the potential energy curves and the various couplings 
that may occur in the A 1 !]^, c 3 II u , and a 3 £+ manifold 
of electronic states. Also, the assignment of the photoas- 
sociation spectrum requires a detailed knowledge of the 
rovibrational levels close to the dissociation limit. At 
present, the experimental data on the A state recorded 
to date 6 - is far from complete. Therefore, in the present 
paper we report a theoretical study of the spectroscopy 
of the strontium dimer in the A X S]+, c 3 II u , and a 3 E+ 
manifold by state-of-the-art ab initio methods. The plan 
of this paper is as follows. In sec. [IT] we describe the ab 
initio electronic structure and quantum dynamical cal- 
culations. We present the numerical results in sec. IIIII 
and discuss at length the accuracy of the present results, 
compare with the available experimental data, and re- 
port predictions for the on-going experiments 3 - 3 -. Finally, 



in sec. IIVI we conclude our paper. 



II. COMPUTATIONAL DETAILS 

A. Ab initio electronic structure calculations 

In the present study we adopt the computational 
scheme successfully applied to the ground and ex- 
cited states of the calcium dimer— ~—, magnesium 
dime r 39 ' 40 , (BaRb) + molecular ionii, and SrYb het- 
eronuclear molecule^ 2 -. The potential energy curves for 
the lowest singlet and triplet excited ungerade states of 
the Sr2 molecule corresponding to the 1 S + 3 P, 1 S + 3 D, 
and 1 S + 1 D dissociation limits have been obtained by a 
supermolecule method: 



V 



2S+1 \Mu(n\ pSM F SM F SM 



(1) 



where Ej^q denotes the energy of the dimer computed 
using the supermolecule method (SM), and E^ 1 , X=A 
or B, is the energy of the atom X. Here, the molecular 
electronic term is denoted by 2S+1 |A| U where S is the 
total electronic spin quantum number, and A the projec- 
tion of the electronic orbital angular momentum on the 
molecular axis. The excited states were calculated em- 
ploying the linear response theory (equation of motion) 
within the coupled-cluster singles, doubles, and linear 
triples (LRCC3) frameworki^i 5 .. Note that in Refs^^ 
the full configuration interaction correction (FCI) for the 
four-electron valence-valence correlation was added on 
top of the linear response result. However, the results 
for Ca ^ 34 ' 35 and more recently on Mg2 show that for the 
states of interest this correction is so small that it can 
safely be neglected. We assume that this is also the case 
for Sr2. 

Transitions from the ground X 1 E+ state to the ^^ 
and 1 H U states are electric dipole allowed. The corre- 
sponding transition dipole moments are given by the fol- 
lowing expression 46 : 



IM (n^X) = (X 1 Y,+ \r i \(n) 1 \A\ u ) 



(2) 



where n numbers the consecutive nonrelativistic disso- 
ciation limits of the 1 |A| U states, and r^, i = x, y or z, 
denotes the ith component of the position vector. Note 
that in Eq. (|2J|, i = x or y corresponds to transitions to 
1 W. U states, while i = z corresponds to transitions to *E+ 
states. In the present calculations the electric transition 
dipole moments were computed as the first residue of 
the coupled-cluster linear response function restricted to 
single and double excitations (LRCCSD) with two elec- 
tric dipole operators 4 ^. Note that in principle the more 
advanced LRCC3 method could be used to calculate the 
transition moments. However, since the intensities in the 
spectra cannot be measured to a high precision, we have 
employed a less accurate method limited to single and 
double excitations. Comparison of the computed and 
measured atomic lifetimes will show below that such an 



approximation is sufficient for the purpose of the present 
study. 

The rovibrational energy levels of the electronically ex- 
cited states of Sr2 are expected to show some perturba- 
tions due to the nonadiabatic coupling between the elec- 
tronic states. Analysis of the potential energy curves, 
cf. sec. IIII1 reveals angular couplings of the a 3 E+ and 
b 3 E+ states with the c 3 II„ state. Therefore, in this work 
we have computed the most important angular coupling 
matrix elements defined by the expression: 



state a 37 ' 47 , and the corresponding spin-orbit matrix ele- 
ments read: 

<pi{R) = (c 3 n u (£ = 0, A = ±l)|i? SO |a 3 S+(S = ±1, A = 0)) , 

(7) 

tp 2 (R) = (c 3 n u (£ = 0, A= ±l)|i? SO |b 3 S+(E = ±1, A= 0)> , 

(8) 

ip 3 (R) = (c 3 n u (£ = 0, A = ±l)\H so \B' 1 U u ) , (9) 



L(n o ri) = ((n) 2S+1 |A| n |L ± |(n') 25+1 |A| n ), (3) 

with L± the ladder operator of the electronic angular mo- 
mentum and n ■<-> n' denoting the coupling between elec- 
tronic states n and n' (here n will stand for the a 3 E+ or 
b 3 E+ states, n' for the c 3 H u state). Note that the elec- 
tronic angular momentum operator couples states with 
the projection of the electronic orbital angular momen- 
tum on the molecular axis A differing by one. In the 
present calculations the angular coupling between the 
triplet states was computed directly from the multiref- 
erence configuration interaction wave functions limited 
to single and double excitations (MRCI). 

Strontium is a heavy atom, and the electronic states of 
the Sr2 molecule are expected to be strongly mixed by the 
spin-orbit interaction. Therefore, the spin-orbit coupling 
and its dependence on the internuclear distance R must 
be taken into account in our analysis of the spectra in the 
mixed singlet /triplet A 1 !]^", c 3 IT u , and a 3 E+ manifold of 
electronic states. We have evaluated the spin-orbit cou- 
pling matrix elements for the lowest dimer states that 
couple to the 0+ and l u states of Sr2 within the MRCI 
framework. In our case the nonrelativistic states that are 
coupled through the spin-orbit interaction to the 0+ sym- 
metry are c 3 IT u , A : E+ and B^+^LII. Other electronic 
states of Sr2 can be omitted from the present analysis 
due to their very weak couplings with the A+c+a mani- 
fold and significant energetic gaps as compared to these 
electronic states. Thus, the most important spin-orbit 
coupling matrix elements for states of the 0+ symmetry 
are given by 

A(R) = 



(c 3 II u (E = ±1, A = Tl)|#so|c 3 n u (E = ±1, A = T l)> 



£i(fl) = (c 3 n u (E = ±1,A= Tl^solA^ 



£ 2 (i?) = (c 3 n u (E = ±1, A = Tl^solB^ 



(5) 



(6) 



where -ffso is the spin-orbit Hamiltonian in the Brcit- 
Pauli approximation 4 ^ and E denotes the projection of 
the electron spin angular momentum on the molecular 
axis. For the l u symmetry the most important cou- 
plings occur between the c 3 II u , a 3 E+, b 3 E+ and B' 1 !!^ 



Ci(i?) = (a 3 E+(E = ±1, A = 0)|iJ so |B' 1 n„) , (10) 

C 2 (R) = (b 3 E+(E = ±1, A = OMsolB' 1 ^) . (11) 

With the spin-orbit coupling matrix elements at hand, 
we can build up the matrices that will generate the po- 
tential energies of the spin-orbit states that couple to 0+ 
and l u symmetry. The matrices for the 0+ and l u states 
are given by: 



yc 



V" 



and 



*(R)-A(R) 

6(A) 
6(A) 



V Al ^(R) 




6 OR) 



v Bl *i(R) 



(12) 



v u 



( v^i 



MR) 
V Ci(R) 



o 

v h3 ^ 
MR) 



MR) 

MR) 

v c3 MR) 

MR) 



Ci(R) 
UR) 
MR) 



v B>1 MR)J 



(13) 
respectively. Diagonalization of these matrices yields the 
spin-orbit coupled potential energy curves for the 0+ and 
l u states. Note that all the potentials in the matrices 
(TT21 and flT3"|) are taken from the LRCC3 calculations. 
Only the diagonal and nondiagonal spin-orbit coupling 
matrix elements were obtained with the MRCI method. 
Once the eigenvectors of these matrices are available, one 
can easily get the electric dipole transition moments and 
the nonadiabatic coupling matrix elements between the 
relativistic states. It is worth noting that here, unlike 
in the case of Ca2, the B X E+ and B^II^ states are in- 
cluded in the model. This is due to the fact that for 
(4) Sr2 the long-range spin-orbit interactions of the c 3 II u 
state with the B and B' states have some significance 
since they are responsible for the existence of very weakly 
bound states located just below the 1 So + 3 Pi threshold 
that were observed in the photoassociation experiment 
by Zelevinsky et al£±. The long-range character of these 
couplings makes the photoassociation of the ultracold 
strontium atoms possible, yielding a non- negligible reso- 
nant 6C^ es R~ 3 interaction for the 0+ and l u potentials at 
the 1 So + 3 Pi asymptote, and a non- negligible relativis- 
tic transition moment from the X 1 E+ ground state. By 

contrast, the coupling with the A' II„ state which was re- 
ported in Ref. 6 was neglected in the present calculations 



TABLE I. Parameters of the analytical fits of the a 3 E+, b 3 E+, c 3 n„, and A X E^ potentials (in atomic units) for Sr2. Numbers 
in parentheses denote the power of 10. 



Parameter 


a E u 


b 3 E+ 


c J n„ 


A^+ 


A 


6.763580846(2) 


3.812370308(5) 


5.806031460(6) 


1.0713571607(3) 


At 


-2.7286531831(2) 


-2.7785352110(5) 


-3.433599894(6) 


-3.1369577995(2) 


A 2 


4.4648343035(1) 


7.599745922(4) 


7.677597973(5) 


4.259633532(1) 


A, 


-3.4042105677 


-9.26259308(3) 


-7.75564134(4) 


-2.828788402 


A 4 


0.1048331384 


4.2607133921(2) 


3.0330564923(3) 


0.08207453705 


a 


0.03722198867 


0.7652092533 


2.102245892 


0.2810137665 


P 


3.14301181557 


0.7516213273 


1.03238202 


1.195331858 


7 


0.07613828056 


0.148486204 


1.349678173(-3) 


0.0362537402 


C12 


-5.318418476(9) 


1.1998965806(11) 


-1.06415514(10) 


7.278665842(11) 


c 5 








-8.649(2) 


c 6 


4.488(3) 


6.750(3) 


3.951(3) 


2.72(3) 


C 8 


1.426(6) 


2.044(4) 


3.521(5) 


2.285(5) 


Cio 


2.321(8) 


1.010(6) 


3.296(7) 


9.223(7) 



since it is asymptotically zero for the l u states of inter- 
est and influences the R~ 6 and higher asymptotics of the 
(l)l u state, as opposed to the SC^ CS R~ 3 asymptotics due 
to spin-orbit coupling with the B' U u state, cf. Ref42 for 
a simple atomic model and Ref. 5 - for a rigorous explana- 
tion. Note finally that the spin-orbit coupling between 
the c 3 II u , a 3 E+, b 3 E+ and B' 1 !^ states also leads to 
states of 0~ and 2„ symmetry. In the absence of strong 
nonadiabatic effects these states are not optically active 
and we do not discuss them here. 



mimic the scalar relativistic effects, 
were described by the ECP28MDF 



In order to 
some electrons 

pseudopotential 5 -^ from the Stuttgart library. Thus, in 
the present study the Sr2 molecule was treated as a sys- 
tem of effectively 20 electrons. In all calculations the 
[8s8p5d4:flg} basis set suggested in Ref^i was used, aug- 
mented by a set of [lslpl<il/3.g] diffuse functions. In 
the calculations of the potentials this basis was supple- 
mented by a set of spdfg bond functions 5 ^ . The full basis 
of the dimer was employed in the supermolecule calcula- 
tions and the Boys and Bernardi scheme was utilized to 
correct for the basis-set superposition error 53 . Ab initio 
calculations were performed for a set of 22 interatomic 
distances ranging from R — 5 to 50 bohr. In the calcula- 
tions of the potential energy curves and transition dipole 
moments the DALTON code 5 -^ was used. All MRCI cal- 
culations were done with the MOLPRO codeS. We would 
like to emphasize that almost all ab initio results were 
obtained with the most advanced size-consistent meth- 
ods of quantum chemistry, LRCC3 and LRCCSD. Only 
the spin-orbit coupling matrix elements and nonadiabatic 
matrix elements were obtained with the MRCI method 
which is not size consistent. Fortunately, all of the cou- 
plings are important in the region of the curve crossings 
or avoided crossings at short interatomic distances, so 
the effect of the size-inconsistency of MRCI on our re- 
sults should not be dramatic. 



B. Analytical fits 

The computed points of the potential energy curves 

were fitted to the following analytical expression: 

V S+1 ^(R) = e-^« 2 J2A t R* + %h{P,R) 



4=0 






n— 3 



(14) 



where {Ai}j =0 , a, (3, 7, and C12 were adjusted to the 
computed points. The damping function f n (/3,R) was 
employed in the form proposed by Tang and Toennies 56 . 
The long-range coefficients {C2n} 5 =3 were not fitted, but 
fixed at the ab initio values taken from Reb^ 3 -. In our 
case the leading long-range coefficient CJj es describes the 
first-order resonant interaction 5 -^ between the Sr( 1 S) and 
Sr( 1 D) atoms, and was also fixed at the ab initio valued. 
For the triplet states the resonant interaction term van- 
ishes identically in the nonrelativistic approximation, so 
C T k os is equal to zero. The parameters of the analytical 
fits of the potentials are reported in Table |TJ 

The transition moments fj,Q^ x (R), spin-orbit cou- 
pling matrix elements, A(R), {&(i?)}? =1 , {d(R)}i=i,2, 
{(^(i?)} 3 =1 , and angular coupling matrix elements, 
L(a o c) and L(b ■H c), were fitted to the following 
generic expression: 



X(R) = X 



'CO I 


(A X 


+ A x R+A x R 2 )e 


-< 


B x e 


n— no 


— f (8 X 


R) 



af R—y x R 2 



(15) 



where X stands for /i"^ x (i?), A(R), {&(R)} 2 =1 , 
{Ci(-R)}f=i, W(R)}U, i(aoc), and I(b«c). The 
leading power in the inverse power expansion of Eq. (|15[) 
depends on the asymptotic multipole expansion of the 
wave functions in the polarization approximation 5 -^— 



and of the appropriate operator and varies between 3 
and 6 for the different quantities X. The atomic values, 
X°°, were fixed as follows: 



TABLE II. Parameters of the analytical fits of the spin-orbit 
coupling matrix elements (in atomic units) for Sr2. Numbers 
in parentheses denote the power of 10. 



A°° = u>?° = 193.68 cm" 1 , 

r , ( 16 ) 

v f = -(°° = ff%/2 = -153.02 cm" 1 , 



too ,.00 /-OO r\ 



(17) 



ACTA^X)=0, Mo°(B-<-X) = V^-3.07a.u. ) (18) 



L°°(aoc) = v / 2, £°°(boc)=0, (19) 



Parameter A(R) 


6 OR) 




&(R) 


X°° 


8.824756(-4) 






-9.860287(-4) 


B x 




0.1281809 




-8.7008164(-3) 


at 




0.255215595 




0.40876725 


At 


0.254210634 


-8.18616724(- 


-2) 


2.96378955(-6) 


At 


-3.332647733(-2) 


9.17710750(- 


-3) 


-3.39493117(-7) 


At 




-4.77060945(- 


-4) 


1.0586325(-8) 


„x 

a 2 


0.7905806772 






-0.92238515 


1 X 




1. 40267678 (- 


-2) 


2.8575636(-2) 


P X 


1.3932739215 


6.8008294 




2.22429216 


x 3 








2.64888569(-2) 


x 4 




-1.684296847 






x 6 


-35.590169218 









and the remaining parameters were adjusted to the ab 
initio points. The value for A°° was derived from the 
experimental positions of the states in the 3 P multiplet 
assuming pure LS coupling, while the values of £2° and 
jUo°(B 4— X) are based on the present ab initio atomic 
calculations. The parameters of the analytical fits of 
the most important spin-orbit coupling matrix elements 
A(R), £i(R), and &(R) are reported in Table HU All pa- 
rameters for other fitted quantities can be obtained from 
the authors upon request. 

Note that fixing our fits at their proper asymptotic val- 
ues is crucial for a proper description of the rovibrational 
transitions near the dissociation threshold. This is in a 
sharp contrast with some potentials fitted to the experi- 
mental data that may not be sufficiently sensitive to the 
long-range tail of the potential and spin-orbit couplings. 



C. Quantum-dynamical calculations 

In the present paper we consider the homonuclear 
bosonic 88 Sr2 molecule. The rovibrational energy levels 
and wave functions for the ground X 1 Ei" state were ob- 
tained by diagonalizing the Hamiltonian for the nuclear 
motion in the Born-Oppenheimer approximation with 
the variable step-size Fourier grid representation^ - — . 
For the ground state, an accurate potential fitted to the 
experimental high-resolution Fourier transform spectra is 
available^ and is used in our calculations. 

Due to the spin-statistical weights for the bosonic 88 Sr2 
molecule, we can limit ourselves to odd values of the ro- 
tational quantum number J and e parity levels^. Rovi- 
brational energy levels for the excited electronic states in 
the A 1 E+, c 3 II u , and a 3 E+ manifold were obtained by 
diagonalizing the following Hamiltonian: 



(H c d ^-A(R) 


6 


6 


-Cj L 


-Cj L 


Cjs 


\ 


?i 


""diag 














V2C JL 


6 





""diag 











V2C JL 


— Cjl 








rr a ^u 
-"diag 





<pi + Cls 





-Cj L 











-"diag 


f2 + Cls 


C2 


Cjs 








(pi + C LS 


¥2 + Cls 


rrc 3 n„ 
-"diag 


<P3 


V 


V2Cj L 


V2C,j L 


Ci 


(2 


¥>3 


^diag / 



where the diagonal term is given by: 



f>(n) (2S + 1) |A| 
^diag 



1 
2/i 



A 



v 



|AU 



(R) 



J{J + 1) + S(S + 1) + L(L + 1) - ft 2 
2fiR 2 



(20) 



A 2 



(21) 



r 



with pr denoting the radial momentum operator, /1 the reduced mass of the dimer, and fl is the projection of 



the total electronic angular momentum on the molecular 
axis, so ft = A + S. The first three rows of the H matrix 
correspond to the states with |fi| = while the last four 
rows to the states with |fi| = 1. The angular (Coriolis- 
type) couplings are defined as: 



where Hv is the so-called Honl-London factor, 



C Jh (R) = - 



Cjs(R) 



[J(J + l)] 1/2 L(R) 



2[iR 2 



[2J(J+1)] 



1/2 



2[iB? 



C -W-^- 



(22) 



(23) 



(24) 



where L(R) stands for L(n f-> n') defined by Eq. ([3]) 
with n and n' properly chosen, the quantum number L 
appearing in Eq. (|21[) is the electronic angular quantum 
number of the excited state atom, and all other symbols 
appearing in Eq. (J2TJ)) are defined in Eqs. (H])-© and 
((7))- (fTT]) . We refer the reader to Ref. 37 for a rigorous jus- 
tification of the above expressions. Note that the term 
2 „i is not rigorously correct, since it results from the 
so-called adiabatic (diagonal) correction for the nuclear 
motion and the above mentioned value is true only in 
the separated atoms limit. At present, there is no ab ini- 
tio electronic structure code that could provide us with 
(even approximate) values of the angular part of the adi- 
abatic correction, so we keep it at its asymptotic atomic 
value. This approximation should work very well for the 
rovibrational levels near the dissociation threshold. If 
the nonadiabatic angular coupling matrix elements are 
small, we can set the Coriolis coupling constants Cjl, 
Cjs, and Cls equal to zero, and the matrix H becomes 
block diagonal with two blocks corresponding separately 
to the 0j~ and 1„ levels: 



(25) 



/ f>c 3 n u 


A(R) 


6 
frA^i 
-"diag 




6 



-"diag 


\ 
/ 


-"diag 


<Pl 



~h 3 y+ 
n diag 

<P2 


¥>2 

o-c 3 n„ 

-"diag 


Ci 

C2 

<^3 


\ 


V Cl 


(2 


<^3 


-"diag 


/ 



(26) 



The line strength in the spectra, S(v'J' <— v"J"), from 
the rovibrational level \v",J") of the ground electronic 
state to the rovibrational level \v', J') of the A 1 !^, c 3 n„, 
and a 3 £ + manifold is given by: 

S(v'J' <- v"J") = (2J' + i)H,j, 



E E 52(Xu»o(v")K(n' <- X)\ Xn ,j,n.{v')) 



O'=0,±1ct= 



-1 n' 



for J' = J" - 1, 



H, 



J'+i 

2,7' + l 

I for /' — J" 

2J' + 1 iU1 J ~ J 1 

■J' r„r. t' — in 

2.]' 



(28) 



for J' = J" 



1. 



and Xu"o( v ") is the rovibrational wave function of the 
ground state, while Xn'j'ii'W) is an eigenfunction of the 
Hamiltonian, either 11 given by Eq. (|20|) . or H°" and 
H 1 " given by Eqs. (|2"5j) and (|2"6")l . if the helicity decou- 
pling approximation (neglecting of the angular Coriolis 
couplings) is employed. 

We also study the positions and lifetimes of Feshbach 
resonances appearing just above the 1 So + 3 Pi dissoci- 
ation threshold. The Feshbach resonances of interest are 
quasi-bound rovibrational levels of the (2)0+ electronic 
state lying above the dissociation limit of the (1)0+ state, 
i.e., in the continuum of this state. Formally, the Fesh- 
bach resonances can be characterized by complex energies 
of the form E r — (i/2)T, where E r denotes the position of 
the resonance state and Y its width which is directly re- 
lated to the lifetime r by r = H/Y. We have determined 
these complex energies by diagonalizing the Hamiltonian 
for the + states, Eq. (1231) . with an imaginary absorbing 
potential Vcap added to the diagonal Hamiltonian terms 
(EH 64 -—: 



,R < Ra, 
,R a <R<Rn 



Vcap(R) 




(27) 



(29) 
where C = 2.62206 and x = (R- i? a )/(i? max - R a )- The 
parameters E m - m , i? a , and i? max were adjusted to obtain 
stable results with respect to small variations of these 
parameters. Approximate positions of the resonances 
were first determined by using the stabilization method 67 
with respect to the size of the grid. With these positions 
at hand, a set of parameters leading to stable complex 
eigenvalues corresponding to the positions and widths of 
the resonances was easily found, yielding R a = 30 bohr, 
i? max = 200 bohr, and E m { n — 30 cm -1 . The knowledge 
of the approximate positions was particularly useful to 
determine the value of E m - m . 



III. NUMERICAL RESULTS AND DISCUSSION 

A. Ab initio electronic structure data 

Before discussing the potential energy curves, we first 
discuss the atomic excitation energies obtained from the 
LRCC3 calculations and the atomic lifetimes. In Table 
IIIII we present calculated excitation energies in compar- 
ison with fully relativistic atomic calculations of Porsev 
and collaborators^ and experimental data. Our pre- 
dicted position of the nonrelativistic 3 P state is 14570.8 
cm , to be compared with the experimental value of 
14704.9 cm" — deduced from the positions of the states 



TABLE III. Excitation energies (in cm ) for the low-lying 
energy levels of strontium atom. 



TABLE IV. Comparison of the present and most recent the- 
oretical and experimental values of the lifetimes of low-lying 
excited states of strontium atom. 



Excited state 


Present 


Reference 68 


Expcriment 0M 


3 Po 


14187.3 


14241 


14317.5 


3 Pi 


14372.7 


14448 


14504.4 


3 P 2 


14762.6 


14825 


14898.6 


3 D! 


18582.9 


18076 


18159.1 


3 D 2 


18637.2 


18141 


18218.8 


3 D 3 


18726.0 


18254 


18319.3 


X D 2 


20650.3 


19968 


20149.7 


x Pi 


21764.3 


21469 


21698.5 



Excited state 



lifetime 



Reference 



in the 3 P multiplet and the Lande rule. For the 1 Pi 
state we obtain 21764.3 cm" 1 , again in very good agree- 
ment with the experimental value of 21698.5 cm" 1 —. For 
the nonrelativistic 3 D state we obtain 18668.8 cm" 1 , to 
be compared with 18255.2 cm" 1 — deduced from the po- 
sitions of the states in the 3 D multiplet and the Lande 
rule. Finally, our term energy for the 1 D 2 state is 20650.3 
cm" 1 , again in a satisfactory agreement with experiment, 
20149.7 cm" 1 —. The accuracy of the atomic spin-orbit 
couplings can be judged by comparing the computed and 
observed splittings of the energy levels in the 3 P and 
3 D multiplets. For the 3 P multiplet, theoretical split- 
tings between the 3 P2 and 3 Pi, 3 P2 and 3 Po, and 3 Pi 
and 3 Po levels amount to 389.9 cm" 1 , 575.3 cm" 1 , and 
185.4 cm" 1 , respectively, to be compared with the ex- 
perimental numbers, 394.2 cm" 1 , 581.0 cm" 1 , and 186.8 
cm" 1 , respectively. A somewhat less good agreement 
is observed for the 3 D multiplet. The theoretical fine 
splittings for the 3 D 3 - 3 D 2 , 3 D 3 - 3 Di, and 3 D 2 - 3 Di states 
read 88.9 cm" 1 , 143.1 cm" 1 , and 54.2 cm" 1 , respectively, 
while the experimental numbers are 100.5 cm" 1 , 160.2 
cm" 1 , and 59.7 cm" 1 , respectively. Finally, we also note 
that the lifetimes of the 3 Pi, multiplet 3 D, 1 D 2 , and 1 Pi 
states of Sr are accurately reproduced. Our calculated 
lifetimes together with the most recent experimental and 
other theoretical results are listed in Table IIVI For the 
1 Pi state we obtained 5.09 ns to be compared with the ex- 
perimental value of 5.22(3) ns2S. For 3 Pi the theoretical 
and experimental numbers are 21.4 /xs and 21.5(2) /is2£, 
respectively. For the D states we observe a slightly worse 
agreement. The theoretical lifetime of the 1 D 2 state is 
0.23 ms, to be compared with the experimental value of 
0.30 ma£i. The same numbers for the average multiplet 
3 D are 2.72 /us and 2.5(2) /j,s~. Such a good agreement 
between theory and experiment for the atoms gives us 
confidence that the molecular results will be of similar 
accuracy, i.e., at worst a few percent off from the exact 
results. 

One of the important issues in ab initio electronic 
structure calculations is the quality of the basis set and 
of the wave functions. To further judge the quality of the 
basis set used in our calculations we have computed the 
leading Cq van der Waals coefficient for the ground X 1 !]^ 



3 Pl 



X D 2 



'■D 



5.09 ns 

5.38 ns 

5.35 ns 

5.22(3) ns 

5.263(4) ns 
21.40 /is 
24.4 fis 
19.0 fis 
21.5(2) fis 

0.23 ms 

0.412(10) ms 

0.30 ms 

2.72 fis 

2.4 fis 

2.5(2) ^s 



Present 
Theory, ReiM 
Theory, Ref.^ 
Experiment, Ref.— 
Experiment, Ref.— 
Present 
Theory, Ref.-^ 
Theory, RetM 
Experiment, Ref.— 
Present 

Experiment, Ref.— 
Experiment, Ref.— 
Present 
Theory, Ref.^ 
Experiment, Ref. 72 



state. This coefficient was obtained by using the explic- 
itly connected representation of the expectation value 
and polarization propagator within the coupled cluster 
metho d 73 ' 74 , and the best approximation XCCSD4 pro- 
posed by Korona and collaborators^. Our ab initio re- 
sult is 3142 a.u. which compares very favorably with the 
value fitted to high-resolution Fourier transform spectra, 
3168(10) a.u. 1 . The agreement between theory and ex- 
periment is better than for most of the other ab initio 
calculations^"—. Comparison between theory and exper- 
iment for the well depth of the ground state X 1 E+ poten- 
tial is somewhat less satisfactory. Our theoretical value 
is 1124.0 cm" 1 , to be compared with the experimental 
result of 1081.64(2) cm"U, i.e., 3.8% too large. How- 
ever, in the case of the ground state interaction, the FCI 
correction for the valence-valence correlation turned out 
to be important. Due to computational limitations, we 
could obtain it only in the Sadlej pVTZ basis^ which 
is comparatively small and does not allow for a better 
accuracy. 

The nonrelativistic potential energy curves relevant for 
the spectroscopy in the A+c+a manifold are plotted in 
Fig. [TJ while the spectroscopic characteristics of these 
states are reported in Table [V] The separated atoms 
energy for each state was set equal to the experimen- 
tal value. Due to the computational limitations in the 
present work we did not consider the B 1 !]^, A n H u and 
B' 1 !!^ states. Fortunately enough, the calculations of 
these potentials were not crucial for our study. Indeed, 
the potential energy curve for the B 1 !]^ state, with the 
correct C^ os R~ 3 asymptotics, could be fitted to the ex- 
perimental data, and is available in the literature*. Ac- 
cording to RefJi, the A n H u state is isolated, and its 
rovibrational levels are only weakly perturbed. Thus, 
it can safely be omitted from our model for the rovibra- 
tional levels of the l u states near the 1 Sq + 3 Pi thresh- 
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FIG. 1. Ab initio potential energy curves for the A E^J 
c 3 IT u , a 3 Si , and b 3 E„ states of the strontium dimer. 




(c 3 n ih — ii- 3 i 




10 15 20 25 30 
R / bohr 



10 15 20 25 30 35 
R/bohr 



FIG. 2. Nonadiabatic matrix elements (left-hand panel) 
and spin-orbit coupling matrix elements (right-hand panel) 
as functions of the interatomic distance R. 



old. However, would this state become interesting from 
an experimental point of view, an accurate potential en- 
ergy curve fitted to the observed spectroscopic transitions 
is available in Refi£, while the important spin-orbit cou- 
pling matrix elements can be obtained from the present 
authors upon request. The second state dissociating into 
the X S + 2 P atoms, the B' 1 !!^ state, was not observed 
experimentally. In our work the potential energy curve 
for this state was approximated with its long-range form 
Cg cs i?~ 3 , with the proper long-range coefficient adapted 
from Ref^ 1 - 3 -. While this is an approximation, it is not a 
crucial one, since the spin-orbit coupling of this state with 
the c 3 II u state is important only at large interatomic dis- 
tances, where it becomes constant, and affects the rovi- 
brational dynamics of the Sr2 molecule only near the in- 
tercombination line 1 So + 3 Pi- 

Let us compare our results with other available ab ini- 
tio data from nonrelativistic calculations 1 —. The spectro- 
scopic constants are listed in Table [V] and compared to 
the results of Ref. 15 . Inspection of Table IV1 shows that 



TABLE V. Spectroscopic characteristics 
relativistic electronic states of Sr2 dimer. 



of the non- 



State 


D e /cm _i 


i? e /bohr 


Ref. 


Dissociation 


a 2-i U 


6298 


7.95 


Present 


i S+ 3 P 




6895 


7.74 


— 






6683 


7.79 


14 




b 3 E+ 


5293 


7.33 


Present 


1 S+ 3 D 




4698 


7.25 


15 






5763 


7.08 


14 




c 3 n„ 


1422 


8.41 


Present 


1 S+ 3 P 




1892 


8.02 


15 






1785 


8.16 


14 




A E u 


8433 


7.54 


Present 


1 S+ 1 D 




5440 


7.12 


15 






9066 


7.28 


14 





the agreement with the data of Czuchaj et alJ^ is not 
satisfactory. For most of the states the computed elec- 
tronic binding energies agree within a few hundred cm -1 
at best, while the differences in the positions of the min- 
ima are 0.4 bohr at worst. The most striking difference 
between the present results and the data of RefJ^ is the 
spectroscopically active A-'-Ej state. Here, the differ- 
ence in the position of the minimum is 0.42 bohr, and 
the difference in D e is as much as 4000 cm -1 . Surpris- 
ingly, the agreement of the present results with the older 
1996 calculations by Aubert-Frecon and collaborators 14 
is comparatively good. Except for the c 3 n„ state, the 
well depths agree within 6% to 8% and the well positions 
within 0.3 bohr at worst. For the c 3 n M state we note a se- 
rious disagreement to all previous results. Our potential 
is considerably shallower and the minimum is shifted to 
larger distances. However, as will be shown in sec. IIIIBI 
the present picture of the interatomic interactions in the 
A+c+a manifold reproduces all features of the available 
experimental data. 

The nonadiabatic and spin-orbit coupling matrix el- 
ements as functions of the interatomic distance R are 
reported in Fig. [5] Note the maximum in the angular 
coupling between the c 3 n„ and the b 3 E+ states. The po- 
sition of this maximum corresponds to the crossing of the 
potential energy curves of these states. Also worth noting 
is the broad maximum of the spin-orbit coupling between 
the A 1 E+ and c 3 II u states. This maximum extends to 
the region where the potential energy curves cross, and is 
responsible for the strong mixing of the singlet and triplet 
rovibrational energy levels. At large distances this par- 
ticular coupling tends to zero, but the coupling between 
the B 1 !]^ and the c 3 n„ states becomes important and 
is responsible for the nonvanishing relativistic dipole mo- 
ment between the ground state and the triplet levels in 
the A+a+c manifold. 

The relativistic 0+ and 1„ potentials as functions of the 
interatomic distance R are depicted in Fig. [3] Their spec- 
troscopic parameters are reported in Table fVll First, we 
note an excellent agreement between the present ab initio 
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TABLE VI. Spectroscopic characteristics of the relativistic 
electronic states of Sr2 dimer. 
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FIG. 3. Relativistic potential energy curves for the lowest 
Oj and 1„ states of the strontium dimer. The dots represent 
the effective empirical potential fitted to the high-resolution 
Fourier transform spectroscopic data of Ref.—. 



potential for the (1)0+ state and the effective empirical 
potential fitted to the high-resolution Fourier transform 
spectroscopic data 6 . The well depths, 2782 cm -1 on the 
theory side and 2790 cm -1 from the fit to the experimen- 
tal data, agree within 8 cm -1 , while the well positions 
agree to within 0.05 bohr. It is gratifying to observe that 
the present ab initio calculations also reproduce satis- 
factorily the position and energy of the avoided crossing. 
Theory predicts the avoided crossing between the (1) and 
(2)0+ potentials at R = 9.1 bohr and V = -1357 cm -1 , 
while the experimental numbers are R — 9.1 bohr and 
V = -1470 cm -1 . We note a substantial disagreement 
between the present ab initio results and those reported 
by KotochigovaiS. The difference in the well depth is 
as large as 946 cm . This means that RefJ^ does not 
predict any interaction between the (1) and (2)0+ states. 
Thus, our result fully confirms the disagreement already 
noticed by Tiemann and collaborators^. For other states 
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FIG. 4. Transition dipole moments between the ground elec- 
tronic state and the excited states of Sr2 in the non-relativistic 
basis (left-hand panel) and the relativistic basis (right-hand 
panel) as functions of the interatomic distance R. 
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16 
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4863 


7.97 


Present 


'So+'Di 




4849 
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16 





the agreement between the two calculations is quite er- 
ratic. For instance, the well depths for the second state 
of 0+ symmetry and the first state of l u symmetry differ 
substantially, while the results for the (2)l u and (3)l u 
states are very close in energy, but the positions of the 
wells are shifted by as much as 0.5 to 0.9 bohr. 

Fig. [4] shows the nonrelativistic electric transition 
dipole moments from the ground state to the 1 £+ states 
and the relativistic transition dipole moments to the 0+ 
states. Inspection of Fig. 0] reveals a broad maximum 
in the transition moment to the A 1 !^ state at distances 
around the minimum of the ground X 1 !]!" state, and de- 
cay to zero at large distances. The transition moment 
to the B^E^ state exhibits a broad minimum in the very 
same region and tends to the atomic value at large dis- 
tances. The two relativistic curves reported in the right- 
hand panel of Fig. @J when superimposed, reproduce the 
transition moment to the A^J state. This is not sur- 
prising, since at small distances the relativistic transition 
moment to the (1)0+ is dominated by the singlet com- 
ponent. Around the crossing between the A-'-Ej and the 
c II U it drops off drastically due to the predominantly 
triplet character of the electronic wave function. Note 
that at large distances it does not decay to zero, but 
rather to a small but constant atomic value reflecting 
the finite lifetime of the atomic 3 Pi state. The oppo- 
site picture holds for the transition moment to the (2)0+ 
state. At small distances this state is dominated by the 
triplet component, and the transition moment is very 
small. Starting from the curve crossing, the (2)0+ state 
is predominantly singlet, and the transition dipole mo- 
ment becomes very large. 



B. Energy levels and rovibrational spectra of Sr_. in the 

A 1 !!*, c 3 n u , and a 3 E+ manifold 

Before comparing our results with the existing exper- 
imental data, let us discuss the validity of the decou- 
pling between the 0+ and l u states. Similarly as for the 



10 



Caz^Z and SrYb^ 2 . molecules, the nonadiabatic angular 
coupling is completely negligible also for the strontium 
dimer, except for a few most weakly bound levels. We 
have calculated the bound states for all possible values 
of the rotational quantum number J by diagonalization 
of the full Hamiltonian matrix (QUI) with the variable 
step Fourier grid representation, and submatrices corre- 
sponding to the 0+ and l u blocks, Eqs. ([251) and (1251) . 
We found that the eigenvalues of the decoupled matri- 
ces differ from the eigenvalues of the full matrix by less 
than 10" 3 cm" 1 , justifying the decoupled representation. 
Actually, experimental data for the rovibrational energy 
levels of the (l)l u state are scarce and limited to three 
levels closest to the x So + 3 Pi dissociation limit, so we 
will focus our discussion on the 0+ levels. 

We have generated all energy levels for the 0+ and l u 
states up to and including J' =219. Comparison of the 
computed energy levels for the 0+ state with the data 
derived from the high-resolution Fourier transform spec- 
troscopy show a root-mean-square-deviation (RMSD) be- 
tween theory and experiment of 10.5 cm -1 . While such 
an agreement between theory and experiment is very 
good for a system with 76 electrons, we decided to adjust 
the ab initio data to the existing experimental data to 
lower the RMSD. This allows us to make reliable predic- 
tions for the ongoing experiments on ultracold strontium 
molecules 3 - 3 -. It turned out that by slightly changing the 
A and C\i parameters by 0.03% and 0.19%, respectively, 
in the analytical expression for the A^+ potential, Eq. 
CH]), we could reduce the RMSD for J' = 1 to 0.64 cm" 1 . 
With the new values of the Aq and C\2 parameters, the 
root-mean-square-deviation of our results for J' < 50, as 
compared to the raw data of Tiemann and collaborators 
(see the supplementary material of Ref. 6 ) is 4.5 cm" 1 . 
We must admit, however, that the present fit is not per- 
fect for very high values of J'. Fortunately, these values 
of the rotational quantum number J' are not of interest 
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FIG. 5. Left-hand panel: Comparison of the original A 1 S^ 
potential with the fitted one to the experimental datar 1 for 
J' = 1. The solid blue line shows the absolute difference, 
AV = Vfittcd-Vao initio, while the green dashed line shows 
the relative difference, ^f = v ^f~ v ^ unite . 100%. Right- 

' V ab initio 

hand panel: the same for the spin-orbit coupling term £2 = 
<c 3 n u |H so |B 1 £+). 



for the ultracold experiments such as photoassociation 
spectroscopy^ or precision measurements leading to the 
determination of the time variation of the electron to pro- 
ton mass ratio 27 ' 28 . To better reproduce the levels near 
the ^o + 3 Pi threshold 2 ^, we have adjusted the spin- 
orbit coupling between the c 3 II u and B 1 E+ states, Eq. 
([6]). Specifically, we varied the parameters X°°, B x , and 
A3 in the fit of &(R)- The adjustments introduced to our 
ab initio data are illustrated in Fig. \5\ Quantitatively, 
the adjusted and ab initio potentials for the A 1 !]^ state 
differ by 290 cm" 1 in the repulsive region, at R = 6.5 
bohr, by -2.6 cm" 1 at the minimum, R = 7.54 bohr, and 
by 0.07 cm" 1 in the long range, at R — 20 bohr. All 
these differences represent at most 2.5% of the original 
ab initio potential. Such an accuracy of the ab initio 
calculations for a system of this size would be purely ac- 
cidental. For the spin-orbit term ^(-R), the adjustment 
results in a change by nearly 50% in the repulsive region 
and less than 0.3% in the asymptotic value. Here, the 
shifts were relatively more important than in the case of 
the A 1 S+ potential, and this can be attributed to the 
lower accuracy of the results from the MRCI method, 
which was employed for the calculation of the spin-orbit 
coupling. It should be stressed, however, that only the 
long-range value of the ^(R) coupling has some signif- 
icance for the described dynamics of the Sr2 molecule, 
as it affects positions of the most weakly bound levels 
below the ^oT^i threshold. Note that the parame- 
ters reported in Tables HI and HT1 are those adjusted to the 
experimental data for J' = 1. 

From Figs. [1] and [3] we expect that some of the rovi- 
brational levels can be assigned to a single state, A-'-E* 
or c 3 n u , while some of them will show a strongly mixed 
singlet/triplet character. This is indeed the case, as is 
nicely illustrated by Fig. El The contour plots of the 
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A 1 !] 



and c 3 II u state populations clearly show that for 



lower values of J' all levels with v' < 18 are predomi- 
nantly singlet levels of the A : E+ state. By contrast, the 
v' = 19 level is the first level that can fully be assigned 
to the c 3 n u state. This is also clear by comparing the 
energies of this particular level in the coupled model and 
in the Born-Oppenheimer approximation: 1363 cm" 1 vs. 
1398 cm" 1 . At higher values of v' the situation becomes 
quite erratic and most of the levels are of strongly mixed 
singlet/triplet character. Only at high values of v' the 
levels can again be assigned, this time to the c 3 II u state, 
although some regions of strong mixing are observed, e.g., 
near v' — 84 or v' — 95. Such a mixing in the rovibra- 
tional levels close to the dissociation limit is very im- 
portant for the ongoing photoassociation experiments of 
ultracold strontium atoms and is discussed in detail in 
Refi 7 .. 

More detailed information on the strongly mixed levels 
is obtained by examing the rotational constants B v > as 
function of the binding energy of the rovibrational level 
\v', J'). This is illustrated in the left-hand panel of Fig. 
[3 showing the first 19 levels from v' = to 18 not to be 
perturbed, with a uniform nearly linear dependence of 
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FIG. 6. Population of the c 3 II u and A 1 S^" components of the 
Oj rovibrational levels lying below the 1 So + 3 Pi asymptote. 
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the rotational constant on the binding energy. Starting 
from v' = 18 the linear dependence is clearly broken, and 
the curve is characterized by irregular behaviour. How- 
ever, it is still possible to distinguish levels which are al- 
most purely triplet. The rotational spacings as functions 
of J'(J' + 1) reported in the right-hand panel of Fig. 
[7] give information about the perturbed levels. Indeed, 
the rotational spacing for v' = shows a purely linear 
character as a function of J' (J' + 1). For v' = 16 this 
linear dependence is broken for high values of J', but the 
nonlinear character is due to the centrifugal distortion, 
and not to the triplet perturbations only. Starting from 
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FIG. 7. Rotational constant of the 0„ rovibrational levels 
for J — 1' (left-hand panel) and rotational spacing, defined 
by AE vJ /(4J - 2) (right-hand panel). 



FIG. 8. Line strengths corresponding to the P (red solid) 
and R (blue dashed) branches in the band i/ — 18 •«— v" = 
27 (upper panel) and v' — 19 •«— v" — 27 (lower panel). 
The spectral lines in the figures are labeled by the rotational 
quantum number J" of the initial state. 



v' = 18, the rovibrational levels can have very different 
character (purely triplet or singlet, strongly mixed), de- 
pending on the value of J'. This is clearly illustrated by 
the sudden changes in the slope of the linear dependence 
of the rotational spacing on J' (J' + 1). 

Figure [8] reports the line strengths for the transitions 
v' = 18 <- v" = 27 and v' = 19 <- v" = 27 for different 
values of the rotational quantum number J" of the ini- 
tial state. The level v" = 27 was chosen since it plays a 
central role in proposed experiments to determine a time 
variation of the electron to proton mass ratio 27 i 28 . In- 
spection of Fig. [5] reveals that the spectral line strengths 



for the 



= 18 



= 27 transitions show a typical 



dependence on the transition frequency, with a slow in- 
crease as a function of J" ' , and a relatively steep and fast 
decrease for large transition frequencies. By contrast, 
the line strengths for the v' = 19 <— v" — 27 transitions 
show a completely different dependence on the transition 
frequency. First we observe a very slow increase of the 
line strength with the increasing J", see the insert and 
note the difference in the intensity scale, and a very slow 
decrease at large values of the transition frequency. Such 
a behaviour is again a signature of strong perturbations 
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from the c 3 n u state that should be observable in exper- 
iment. 



C. Feshbach resonances above the 1 So + 3 Pi asymptote 

The rovibrational states that could be assigned to the 
(2)0+ state lying above the 1 So + 3 Pi dissociation limit 
are embedded in the continuum of the (1)0+ state, and 
thus are Feshbach resonances, i.e., quasi-bound states 
with finite lifetimes. Two decay mechanisms are pos- 
sible depending on the strength of the coupling with the 
continuum. The resonances can decay to the ^o + 3 Pi 
continuum, i.e., dissociate into atoms (predissociation), 
or, if the lifetime is long enough, decay to bound or con- 
tinuum states of the electronic ground state (spontaneous 
emission). Our computed positions, widths, and lifetimes 
of a few low lying Feshbach resonances for J' — I are re- 
ported in Table IVIII The width of the resonances varies 
quite strongly, between almost 1 cm^ 1 for the very broad 
one to «0.001 cm -1 for the very narrow one. The cor- 
responding lifetimes also vary considerably, from the pi- 
cosecond to the nanosecond scale. This implies that some 
of the resonances should be observable in high-resolution 
Fourier transform spectroscopic experiments. 



In Table IVIII we have numbered the resonant states 
by the vibrational quantum number v' corresponding to 
bound vibrational states above the 1 So + 3 Pi dissocia- 
tion limit that result from the diagonalization of the H°" 
Hamiltonian. For </' = 1 we have found 110 bound lev- 
els of the + symmetry which lie below the 1 So + 3 Pi 
asymptote, assigning them vibrational quantum numbers 
ranging v' = to v' = 109. Thus the first quasi-bound 
level reported in Table IVTll located above the ^o + 3 Pi 
threshold has the vibrational quantum number v' = 110. 
Only eigenvalues that were stable with respect to the grid 
size were selected, and the exact positions and lifetimes 
were determined with the complex absorbing potential 
of Eq. (f2"Tj|). It turns out that the width/lifetime is a 
very sensitive function of the strength of the spin-orbit 
coupling. To illustrate this point, we report the widths 
of the resonant states as a function of the triplet state 
population of the rovibrational wave function in the left- 
hand panel of Fig. [S] The narrow resonances are found 
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FIG. 9. Widths of the Feshbach resonances as a function of 
the population of the c 3 II u state and scaling parameter A of 
the spin-orbit term (c 3 il„|Hso|A 1 E + ). 



to be dominated by the singlet component of the wave 
function, for which the coupling through the spin-orbit 
interaction is weak. By contrast, those resonances that 
are almost purely triplet states are very broad with short, 
picosecond lifetimes. This is easily rationalized in terms 
of a state with a predominantly singlet character decay- 
ing more slowly to the triplet continuum than an (almost) 
purely triplet state. Visibly, the spin-orbit coupling is 
responsible for the predissociation decay of the Feshbach 
resonances, i.e., a smaller spin-orbit interaction favours 
longer- lasting quasi-bound states, cf. v' — 114 reported 
in Table EIIl 

It was shown in Ref. 7 - that the widths of the reso- 
nances are very sensitive to the quality of the ab ini- 
tio data. Since the existence of low-lying Feshbach res- 
onances is mostly due to the spin-orbit coupling £i(i?), 
between the c 3 II K and A 1 S + states, we have checked how 
scaling this particular term by a parameter A ranging 
from 0.9 to 1.1, A • £,i(R), affects the results. This is il- 
lustrated in the right-hand panel of Fig. [J3 Inspection 
of this figure shows the width to be indeed a rather sen- 
sitive function of the coupling, and a 10% change in the 
coupling may result in a change of the width by a fac- 
tor of two. However, since scaling within ±10% does not 
change the order of magnitude of the lifetimes, our con- 
clusions concerning a possible experimental observation 
of the Feshbach resonances in the A 1 ^, c 3 II u , and a 3 S + 
manifold remain valid. 



TABLE VII. Positions (relative to the 1 S + 3 Pi threshold), 
widths, and lifetimes of the Feshbach resonances of 88 Sr2 
above the So + 3 Pi dissociation limit. 



IV. SUMMARY AND CONCLUSIONS 



v' 


Er/cm- 1 


r/cm" 1 


r/ns 


110 


49.8 


0.8761 


0.0037 


111 


112.2 


0.4658 


0.0067 


112 


174.4 


0.1924 


0.0162 


113 


236.3 


0.0410 


0.0759 


114 


297.9 


0.0008 


3.8821 


115 


359.2 


0.0638 


0.0488 



We have reported a theoretical study of interatomic 
interactions and spectroscopy of the strontium dimer in 
the A 1 !]^, c 3 II u , and a 3 S + manifold of electronic states. 
The spectroscopic characteristics of the potential energy 
curves for the nonrelativistic A 1 E + , c 3 II u , and a 3 E+ 
states of Sr 2 and of the relativistic states of + and l u 
symmetry deviates significantly from most of the previ- 
ous ab initio result o 15 ' 16 . This is particularly true for the 
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optically active A-'-E^ state and its Hund case (c) analog, 
(1)0+. By contrast, the present spectroscopic character- 
istics of the (1)0+ state are in a very good agreement with 
the experimental values deduced from the high-resolution 
Fourier transform spectroscopic data 6 . The rovibrational 
energy levels corresponding to the spin-orbit coupled 
(1)0+ state dissociating into x So + 3 Pi atoms lying be- 
low the avoided crossing with the (2)0 + state dissociat- 
ing into ^o + 1 T>2 atoms are almost unperturbed, and 
the corresponding energies are very close to the ener- 
gies obtained in the Born-Oppenheimer approximation. 
The rovibrational levels of the (1)0 + state lying above 
the avoided crossing with the (2)0 + state are all heav- 
ily perturbed by the rovibrational states of the c state. 
These perturbations are exclusively due to the spin-orbit 
interaction. In all cases, the nonadiabatic effects due 
to the Coriolis couplings were shown to be negligible, 
with the exception of a few least bound levels. We have 
located several quasi-bound Feshbach resonances lying 
above the 1 So + 3 Pi dissociation limit. Their lifetimes 
suggest that they should be observable in high-resolution 
spectroscopic experiments. 

Overall, our results reproduce (semi) quantitatively the 
experimental data observed thus far—. Our spectro- 
scopic predictions for on-going and future experiments 
concerning the photoassociation of ultracold strontium 
atoms 26-33 and precision measurements of the time varia- 
tion of the electron to proton mass ratio 27 i 28 are reported 
elsewhere 7 - 7 -. 
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